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ABSTRACT: 


Using  laser  frequency-modulation  (I'M)  spectroscopy  and  persistent  spectral  holes, 
time-resolved  phase-sensitive  probing  of  ballistic  phonons  generated  by  an  absorbed 
Nd.YACJ  laser  pulse  is  accomplished  in  the  interior  of  a  solid  sample.  Measurement  of  the 
time  dependence  and  sign  of  a  propagating  stress-strain  lii  illustrated  using  spectral  holes 
in  the  inhomogcncously  broadened  607  nm  color  center  aU  ion  in  x-irradiated  NaE  at 
liquid-helium  temperatures.  By  examining  the  dependence  of  the  observed  phonon 
timc-of-flight  data  on  the  polarization  of  the  probing  light,  the  position  within  the  sample, 
and  the  phase  of  I'M  detection,  an  identification  of  the  acoustic  polarizations  of  the 
propagating  phonons  may  be  made.  The  effects  of  phonon  focusing  and  mode  conversion 
upon  reflection  must  be  taken  into  account  to  complete  the  identification.  Along  with  the 
ability  to  determine  the  sign  of  the  acoustic  disturbance,  this  experiment  features  a  strain 
detection  limit  of  4  x  10  "  at  a  time  resolution  of  50  ns. 


I.  Introduction 


The  study  of  the  propagation  of  nonequilibrium  phonons  in  solids  has  provided  fertile 
ground  for  the  observation  of  several  novel  physical  effects'.  I  sing  heat-pulse  generation  in 
constantan  films  and  detection  with  superconducting  bolometers.  7  a  fascinating 
phonon-focusing  effect  was  observed  due  to  the  intrinsic  elastic  anisotropy  in  cubic  crvstals 

4.5 

Phonon  focusing  has  led  to  the  development  of  a  variety  of  powerful  techniques  for 
imaging  laser-generated  phonon  pulses6,7,  including  the  studv  of  phonon  storage  effects  in 
the  "hot  spot"  near  the  generation  region  s \ 

In  contrast  to  surface  detection  with  ultrasonic  transducers  and  supconducting 
bolometers,  several  optical  methods  for  the  studv  of  phonons  have  appeared  which  allow 
probing  of  phonon  dynamics  and  propagation  within  the  bulk  of  a  solid  material.  In  an  carlv 
experiment10,  the  R-linc  fluorescence  of  Cr’ •  ions  in  ruby  was  used  as  a  monitor  of  29 
cm  1  phonons.  This  experiment  and  others  allowing  detection  of  high-frequency  phonons 
by  phonon-induced  fluorescence  have  been  reviewed  recently".  In  recent  experiments, 
enhanced  optical  dephasing  as  measured  by  optical  free-iiuluction  decay  12  and  modulation 
of  two-pulse  photon-echo  decays  n  have  both  been  shown  to  be  sensitive  detectors  of  higti- 
and  low-  frequency  phonons. 

Persistent  spectral  holes  also  provide  a  sensitive  optical  detector  of  phonon  wave 
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packets  .  In  persistent  spectral  hole-burning,  the  profile  of  an  inhomogcneouslv  broadened 
optical  absorption  of  a  defect  in  a  solid  at  low  temperatures  is  modified  bv  direct 
photochemistry  or  by  photophysical  modification  of  the  local  environment  of  the  absorbing 
defect  .  flic  narrow  spectral  width  of  the  resulting  persistent  spectral  hole  (PSII)  allows 
interaction  with  acoustic  waves  with  quite  high  sensitivity.  In  the  first  experiment  of  this 
type'  ,  an  ultrasonically-gcncratcd  phonon  pulse  was  shown  to  lead  to  strong 
(phase-insensitive)  modulation  of  a  spectral  hole  with  5  u s  time  resolution.  Hv  careful  control 
of  an  impressed  ultrasonic  field,  the  resulting  shifting  and  splitting  of  a  spectral  hole  was 


_  1_ 

shown  to  lead  to  both  phase-insensitive  and  phase-sensitive  '  detection  techniques  for 
spectral  holes  themselves.  In  recent  studies,  phonon-induced  changes  of  a  spectral  hole  and 
anti-Stokes  phonon  sideband  methods  have  been  used  to  study  heat-pulse  propagation  in  thin 

1 X 

polymer  films  ' . 

In  this  paper,  the  powerful  method  of  laser  frequency-modulation  spectroscopy  (I  MS) 
is  used  to  timc-rcsolvc  the  changes  in  a  persistent  spectra!  hole  with  50  ns  time  resolution 

caused  by  the  passage  of  laser-generated  ballistic  phonon  wave  packets  through  the  detection 
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region.  In  FMS"  ,  a  frequency-modulated  laser  beam  is  converted  to  an  .amplitude-modulated 
beam  by  the  narrow  absorption  and  dispersion  due  to  the  spectral  hole.  The  method  is 
sensitive  to  the  phonon  phase,  that  is.  signals  of  different  sign  arc  produced  bv  dilation  and 
compression.  Since  the  detection  operates  at  a  high  rf  frequency  far  above  the  laser 
low-frequency  noise  spectrum,  the  technique  can  operate  near  the  shot-noise  limit. 

To  illustrate  the  technique,  the  607  nm  /ero-phonon  optical  absorption  of  aggregate  color 
centers  in  NaF  was  chosen  as  the  inhomogcneously  broadened  optical  transition  capable  of 
PS  1 1  formation.  Timc-of-flight  data  for  phonon  wave  packets  generated  by  the  absorption 
of  Nd:YAG  laser  pulses  in  a  metal  film  on  the  NaF  crystal  have  a  rich  structure,  requiring 
a  fairly  detailed  analysis.  In  the  next  section,  the  experimental  methods  used  to  record 
laser-generated  phonon  timc-of-flight  signals  me  described.  In  Section  III,  the  theory  of 
phonon-induced  splitting  and  shifting  of  the  optical  transition  is  reviewed.  In  Section  IV, 
experimental  results  on  phonon  detection  and  probe  polarization  dependence  arc  presented. 
The  technique  is  demonstrated  to  be  sensitive  to  strains  as  small  as  4x  10  \  In  Section  V. 
the  dependence  of  the  timc-of-flight  data  on  laser  position  along  with  the  theory  of  See.  Ill 
arc  utilized  to  identify  the  acoustic  polarization  of  the  various  phonon  peaks.  Since  this 
technique  is  able  to  measure  the  stress  and  strain  of  ballistic  phonon  packets  in  the  interior 
of  a  solid  with  high  time  resolution,  it  should  prove  to  be  useful  lor  the  study  of  phonon 
dynamics  and  propagation  in  a  variety  of  materials. 


II.  Experimental 


Figure  I  shows  the  sample  configuration  (or  laser  generation  and  phase-sensitive 
time-of-flight  detection  of  phonons.  The  material  studied  in  these  experiments  is  single  crystal 
N’aF  containing  aggregate  color  centers.  Ballistic  packets  ok  acoustic  phonons  arc  generated 
bv  Nd:YAG  laser  pulses  incidetit  on  a  thin  him  of  (Y  that  has  been  evaporated  on  one  face 
of  the  NaF  crystal.  Phonons  propagating  a  wav  from  the  (Y  film  into  the  NaF  crystal  arc 
detected  inside  a  selected  crystal  volume  by  a  probe  laser  beam.  The  basic  physical 
phenomenon  used  in  the  detection  process  is  that  phonon  stress  and  strain  fields  shift  and 
split  the  frequencies  of  the  color  center  absorption  lines.  For  sensitive  detection  of  the 
absorption  frequency  shifts,  the  probe  laser  is  fust  used  to  burn  a  narrow  persistent  spectral 
hole  (PSII)  in  the  inhomogcncously  broadened  607  nm  color  center  zero  phonon  line  (ZPF) 

is.  -n-2 1  jjlc  pP|0nons  arc  then  detected  as  small  frequency  shifts  of  the  PSII  using  wide 
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bandwidth  laser  frequency-modulation  spectroscopy  (IMS)"  .  Notice  that  both  the 
phonon-generation  and  probe  laser  beams  may  be  displaced  to  examine  the  propagation 
behavior  of  phonons  inside  the  crystal.  T  he  ('artesian  coordinate  system  shown  in  I  ig.  I. 
with  origin  centered  on  the  Cr-coatcd  end  face  and  axes  along  <  100  -•  cubic  crystal 
directions,  will  be  used  throughout  this  paper  to  describe  the  positions  of  the  laser  beams 
and  the  components  of  the  phonon  stress  tensor. 

The  sample  consisted  of  a  single  crystal  of  NaF,  cleaved  and  polished  along  {100}  crystal 
planes  with  dimensions  of  lxl\2  cm.  Color  centers  were  produced  by  exposing  one  of  the 
Ixl  cm  crystal  faces  to  (  u  x-rays  (70  mA  at  70  kV)  for  705  hours  at  room  temperature. 
The  crystal  acquired  a  gradient  of  pinkish  brown  color  from  the  irradiated  end.  A  0.7  /mi 
layer  of  Cr  was  evaporated  on  the  face  opposite  the  irradiated  face.  For  convenience,  the 
NaF  color  center  aggregate  with  ZPF  at  607  nm  was  chosen  as  the  inhomogcncously 
broadened  optical  transition  capable  of  PSII  formation  ’.  I  he  607  nm  ZPI.  had  an  internal 
transmittance  of  00  %  through  I  cm  near  the  middle  of  the  crystal  at  l.o  K,  as  measured 
with  a  tungsten  lamp.  7/4  m  monochromator,  and  (iaAs  photomultiplier  tube. 
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I'igurc  2  shows  the  experimental  apparatus  lor  phase-sensi'iv  e  time-of-llight  detection  of 
ballistic  phonon  packets.  Hie  subsystems  of  the  apparatus  are  a  frequency-modulated  CW 
dye  laser  probe  beam,  a  pulsed  Nd:YAG  laser  beam  used  to  generate  phonons  in  the  sample 
held  at  1.6  K  in  an  optical  immersion  cryostat,  and  high-speed  lascr-frcquencv-modulation 
detection  electronics.  Details  for  each  of  these  subsystems  w  ill  be  presented  first,  followed 
by  a  brief  elaboration  for  certain  components  and  the  procedure  used  to  obtain  phonon 
time-of-flight  signals. 

A  single-frequency  R6('r  dvc  laser  (Df.,  Coherent  596-2!)  beam  with  a  jitter-limited 
linewidth  of  3  Mil/,  passes  through  a  Faraday  optical  isolator  (Ol),  an  amplitude  stabilizer 
(AS),  and  a  FiTaO,  electro-optic  phase  modulator  (FOM.  Crystal  Technology  Model  6100) 
driven  at  250  Mil/.  The  phase  modulator  produces  a  frequency-modulated  laser  beam  which 
is  reflected  from  a  confocal  Fabry  Perot  ctalon  (I  P)  with  7.5  (HI/  free  spectral  range,  is 
attenuated  with  fixed  and  variable  attenuators  (VA),  and  can  be  blocked  with  a  mechanical 
shutter  (SI I)  while  performing  diagnostics  with  the  I  P.  Opening  the  shutter  allows  the  dye 
laser  beam  to  pass  through  a  high-speed  acousto-optic  modulator  shutter  (AOM,  Isomet 
I205-603D)  and  to  be  focused  into  the  sample  (S)  with  a  lens  fl).  The  transmitted  probe 
beam  is  detected  with  a  high  speed  Si  avalanche  photodiode  (API).  Analog  Modules  FNC'A, 
275  Mil/  bandwidth,  -210V  bias).  The  laser  beam  diameter  varies  from  75  to  100  /nn  in  the 
sample. 

A  Q-switchcd  Nd.'YAG  laser  (YAG)  is  operated  at  1.06  /mi  and  produces  8  ns  long  pulses 
at  a  repetition  rate  of  6  II/.  The  YAG  laser  beam  is  split  with  a  glass  beam  splitter  (IFS) 
to  produce  detection-triggering  and  phonon-generation  beams.  The  reflected  S  AG  laser  beam 
is  converted  to  an  electronic  trigger  by  a  300  Mil/  bandwidth  Si  p-i-n  photodiode  (PD). 

The  transmitted  beam  is  expanded  and  apertured  to  <>.l  cm  for  uniform  illumination,  and  is 
attenuated  to  I  - 1 0  /j.l/pulsc  with  a  polarizer  (P)  before  hitting  the  sample.  Since  the  APD 
can  also  detect  YAG  laser  scattered  light,  any  scattered  N  AG  radiation  is  blocked  with 
several  short  pass  interference  filters  (I  ). 


The  AOM  in  the  probe  beam  is  used  as  a  fast  shutter  to  expose  the  sample  to  the 
probing  light  for  only  a  short  time  about  the  YAG  firing  time  (the  reason  for  this  is  explained 
below).  The  AOM  needs  to  be  turned  on  some  time  before  the  YAG  laser  fires.  Since  the 
YAG  laser  lamps  fire  before  the  YAG  reaches  threshold,  the  lamp  signal  is  used  to  trigger 
a  digital  delay  generator  (DPG),  the  output  of  which  is  used  to  open  the  AOM.  .litter  in 
the  lamp-to-lascr  firing  time  is  unimportant,  since  the  transient  digitizer  is  triggered  by  the 
laser  pulse  itself. 

Standard  laser  FMS  techniques  arc  used  to  sensitively  detect  the  time  variation  of  a 
persistent  spectral  hole  in  a  20  MHz  bandwidth  27,2  .  I  lie  sensitivity  results  from  detection 
at  a  frequency  above  the  laser  excess  noise  spectrum  where  the  laser  noise  is  mostly  shot 
noise.  An  RF  source  operating  at  \>R1 .- =  250  MHz  and  4-7  dBm  is  split  with  a  -10  dB 
directional  coupler  (DC),  the  split  signal  of  which  is  phase  shifted  (PS),  amplified  (AI.  +32 
dB),  and  used  to  drive  the  FoM.  The  through  port  of  the  DC  drives  the  local  oscillator  (I. O) 
port  of  the  RF  mixer  (M,  Minicircuits  ZFM-3).  As  is  usual  in  FMS.  the  frequency  spectrum 
at  the  output  of  the  FOM  consists  of  a  carrici  (the  dye  laser  frequency)  and  turn  sidebands 
spaced  above  and  below  the  carrier  by  ±  vR,  with  one  sideband  in  phase  with  the  carrier  and 
the  other  sideband  180°  out  of  phase  with  the  carrier.  When  heterodyne  mixing  of  the 
sidebands  with  the  carrier  occurs  in  the  API),  a  180°  phase  difference  results  between  the 
two  beat  signals  at  vRF.  With  the  sidebands  balanced,  the  beat  signals  cancel  and  there  is 
no  RF  photocurrent  from  the  API)  at  vK|.  If  a  sharp  spectral  feature  such  as  a  persistent 
spectral  hole  with  width  on  the  order  of  vRF  is  probed  by  the  sidebands,  the  resulting 
imbalance  in  the  sidebands  is  translated  into  a  photocurrent  in  the  API)  at  vRF.  The  detected 
RF  signal  is  amplified  by  A2  (  +  34  dB  gain.  2-500  MHz),  and  sent  to  the  RF  (R)  port  of 
the  mixer.  At  the  IF  (I)  port  of  the  mixer,  a  20  Mllz  low  pass  filter  (1  PIT)  blocks  sum 
frequencies,  and  passes  the  envelope  of  the  RF  signal  which  emerges  from  the  API)  in  a  20 
MHz  bandwidth  about  vRF.  The  (time-dependent )  envelope  signal  is  further  amplified  in  A3 
(gain  of  20,  0-150  Mil/)  and  sent  to  either  (I)  a  transient  digitizer  (  I  D)  to  record  the  phonon 
timc-of-flight  signals  with  the  laser  frequency  held  fixed,  or  (2)  a  2  kHz  low  pass  filter  (I.PI  2) 


and  an  oscilloscope  (OS(')  to  obtain  a  I'M  spectrum  of  the  hole  shape  bv  tuning  the  laser 
frequency  about  the  hole. 

Before  phonon  time-of-flight  spectra  arc  obtained,  a  PSII  is  burned  in  the  sample 
absorption  spectrum.  The  probe  laser  is  tuned  to  the  red  of  the  607  nm  /PI.,  the  positions 
of  the  laser  beams  arc  selected  in  the  sample,  and  the  crystal  is  oriented  normal  to  the  probe 
beam.  I  he  probe  laser  polarization  is  selected  with  a  Fresnel  rhomb  and  the  laser  is  tuned 
to  a  position  near  the  center  of  the  /PI.  at  a  wavelength  between  607.145  and  607.155  nm. 
l  hc  laser  vacuum  wavelength  is  measured  with  a  wavemeter.  A  PSII  is  then  burned  into  the 
ZPF  spectrum  with  the  dye  laser  frequence  fixed.  PM  off.  ami  (lie  AOM  on.  The  hole 
burning  time  is  typically  5  to  7  min  using  2fH)  »\V  in  a  beam  diameter  of  75  /mi. 

After  a  hole  is  burned,  the  PSII  PM  spectrum  is  displayed  on  the  DSC  in  order  to  adjust 
the  RF  phase  and  select  a  laser  frequency  for  phonon  detection,  l  hc  adjustment  process 
requires  about  a  minute,  and  so  the  probe  laser  intensity  is  reduced  to  7  /t\V  to  reduce 
subsequent  hole-burning  during  the  set-up  process.  At  this  power  level,  the  2  kHz  low  pass 
filter  bandwidth  of  FPF2  is  adequate  for  recording  the  PSII  IM  spectrum  with  a 
signal-to-noisc  ratio  (SNR)  of  about  10  in  one  0. 2.5  s  scan. 

Figure  3  shows  PSII  shapes  in  the  Nal  607  nm  /PI  at  1.6  K  as  measured  with  IMS 
by  scanning  the  dye  laser.  Figure  5  (a)  is  the  absorptive  or  S,  phase  (in  phase  with  the  RF 
source)  and  (b)  is  the  dispersive  or  S2  phase  (quadrature  phase)  PSII  frequency  modulated 
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spectrum  .  lhc  S|  phase  spectrum  in  Pig.  3  (a)  contains  two  copies  of  the  absorption 
hole  shape,  each  produced  a s  an  I'M  sideband  is  tuned  through  the  PSII.  Modeling  the  PM 
line  shape  reveals  that  these  holes  have  hole  widths  (IAVIIM)  on  the  order  of  the  modulation 
frequency  of  250  Mil/..  I  his  hole  width  is  somewhat  larger  than  the  true  low-intensity  width 
of  45 ±6  MIlz  due  to  the  need  to  use  a  iclativeh  tight  locus  in  the  sample.  The  S_,  phase 
spectrum  is  composed  of  three  overlapping  copies  of  the  PSII  dispersion,  one  centered  at  the 
hole  burning  frequency,  and  two  others  displaced  by  with  half  the  amplitude  of  the 

central  feature.  The  depth  of  the  absorption  hole  producing  the  PM  lineshapes  in  Fig.  3  is 


calibrated  using  the  si/c  of  the  reflection  dip  of  the  eialon  and  the  etalon  IM  spectrum. 
Holes  burned  near  the  middle  of  the  sample  tvpicalh  i  (insist  of  about  a  5'\,  change  in 
transmission. 

To  obtain  plionon  timc-of-flight  signals,  the  Rl  phase  of  detection  is  selected  to  be  either 
Si  or  S?,  the  etalon  signals  arc  tuned  away  from  the  llS 1 1 .  and  the  laser  rre«|uency  is  fixed. 
In  Sj  phase,  the  laser  frequency  is  fixed  at  the  high  slope  recion  at  the  center  of  the  hole 
spectrum  (f  ig.  3  (b)|.  In  S,  phttsc,  the  laser  frequency  is  fixed  at  the  negative  peak  of  the 
negative  copy  of  the  absorption  hole  (lag.  3  (an.  With  no  modulation  of  the  hole  shape, 
the  output  of  the  mixer  remains  constant  at  zero  (for  Sn  or  at  the  negative  peak  value  (for 
S,).  When  phonon  pulses  propagate  thr<  ugh  the  hole-burned  volume,  the  PSII  splits  under 
the  influence  of  phonon  stress  and  strain,  and  the  hole  modulation  signals  arc  recorded  and 
averaged  on  the  I'D. 

As  is  well-known  for  spectral  holes  in  color  center  /PI  s,‘ '  excessive  reading  intensity 
can  cause  distortion  of  the  PSII  shape  by  additional  lude-burning.  At  the  same  time,  since 
the  fast  phonon  timc-of-flight  signals  are  recorded  with  a  20  Mil/  detection  bandwidth, 
higher  reading  powers  arc  required  to  maintain  SNR  compared  to  the  case  with  2  kll/ 
bandwidth.  To  satisfy  these  conflicting  requirements,  the  acousto-optic  modulator  (AOM) 
is  used  as  a  shutter  to  transmit  the  20  /AV  d\e  laser  beam  to  the  sample  only  dining  the 
phonon  timc-of-flight  detection  period.  This  icduces  the  duty  cycle  to  (loo  /<s)(6  II/)  = 
0.0016,  where  100  /ts  is  the  longest  recording  interval  used,  latch  timc-of-flight  signal  is 
typically  averaged  over  1500  YA(i  laser  shots,  for  a  signal  integration  time  less  than  0.15 
seconds  per  averaged  trace,  but  the  total  experiment  time  is  nearly  *1  minutes.  In  this 
manner,  a  higher  phonon  timc-of-flight  SYR  is  achieved  with  lower  total  probe  laser  flux 
and  insignificant  additional  hole  burning. 


III.  Detection  theory 


Phonon  propagation  experiments  frequentlv  employ  pulsed  laser  heating  of  a  metal 
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surface  to  generate  phonon  pulses  "  In  the  nonahlative.  low-poucr  or  thcrmoclastic 

regime  the  laser  radiation  is  absorbed  within  a  lew  hundred  A  of  the  metal  surface.  The 
microscopic  view  at  low  temperatures  is  that  the  absorbed  laser  energy  quicklv  converts  to 
lattice  vibrations  or  phonons  with  a  uniform  distribution  of  phase  velocity  directions.  For 
the  YACi  pulse  llucnccs  used  in  these  experiments,  the  surface  temperature  is  estimated  to 
be  higher  than  15  K.  which  results  in  a  distribution  of  phonon  frcc|uencics  up  to 
approximately  0.5  I'll/.  This  frequency  is  well  below  the  optic  phonon  spectrum  and  in  the 
range  of  acoustic  phonons  for  Nab.  The  laser-pulse-generated  acoustic  phonons  propagate 
into  the  crystal  in  three  separate  expanding  shells  or  wave  packets,  one  composed  of 
longitudinal  and  the  other  two  composed  of  transverse  acoustic  phonons  Since  the 

phonons  contributing  to  the  ballistic  phonon  pulses  we  observe  arc  of  fairly  low  frequency 
and  in  the  linear  dispersion  region,  we  may  in  a  reasonable  approximation  focus  on  only  one 
type  of  low-frcqucncv  phonon  at  a  time.  This  is  a  result  of  the  fact  that  till  the  phonons 
of  interest  in  a  particular  polarization  propagating  in  a  fixed  direction  have  the  same  velocity 
characteristics,  flic  resulting  total  signal  will  be  a  superposition  of  the  physical  effects  from 
all  the  phonons  present. 

The  basic  idea  behind  detecting  a  phonon  with  a  persistent  spectral  hole  is  that  the 
phonon  is  a  travelling  stress-strain  field  and  that  the  spectral  hole  splits  and  shifts  under  the 
influence  of  stress  and  strain.  The  .ability  of  the  probe  laser  beam  to  phase-sensitively  detect 
the  perturbations  tine  to  the  local  stress-strain  field  depends  upon  the  detection  bandwidth 
as  well  as  the  relative  si/e  of  the  probe  beam  diameter  and  the  phonon  wavelength.  In  these 
experiments,  the  probe  beam  diameter  is  smaller  than  one-half  the  minimum  acoustic 
wavelength  expected  for  20  Mil/  phonons  in  NaF;  thetefore.  the  ability  to  follow  the  time 
dependence  (and  sign)  of  phonons  up  to  the  bandwidth  limit  of  20  Mil/  is  expected. 


In  the  standard  anisotropic  continuum  limit,  a  wave  equation  lor  the  displacement  of 
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volume  elements  in  the  solid  has  plane  wave  solutions  of  the  form 

u(r .  t)  —  u„  cxp[  i(k  •  r  --  ot)l.  (I) 

where  ti„  is  the  polarization  of  the  acoustic  wave,  k  is  its  wave  vector,  and  o  is  its  frequency. 
There  arc  three  difTcrcnt  phase  velocity  surfaces  vp  —  (o/k)k  in  k  space  where  k  is  a  unit 
vector  in  the  direction  of  k.  The  three  surfaces  correspond  to  difTcrcnt  orthogonal 
displacement  polarizations.  The  phase  velocity  surface  with  polarization  u~  ncarlv  parallel 
to  vr  is  known  as  the  quasi-longitudinal  surface  and  the  other  two  a  re  the  quasi-shear  or 
quasi-transverse  surfaces.  The  strain  tensor  for  the  displacement  wave  in  Tq.  I  is  formed 
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from  the  symmetric  part  of  the  gradient  tensor  '  ‘  and  has  (real)  components 

r,j  =  — y  (tijkj  +  ttjk,)  sinfk  •  r  -  at).  (2) 

Hquivalcntlv,  the  stress  in  the  wave  may  he  computed  through  the  elastic  stiffness  tensor  as 

=  *  mw'Ti*  0) 

where  summation  over  repeated  indices  is  assumed.  It  is  useful  to  work  with  the  stress  tensor 
since  a  connection  can  directly  be  made  to  previous  uniaxial  stress  calculations  11  and 
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experiments  for  the  607  nm  NaT  color  center. 

following  Kaplyanskii,13  the  absorption  frequency  shift  for  a  single  defect  under 
externally  applied  stress  may  be  written  in  a  linear  approximation  as 

A  =  Axx^xx  f  Ayy^vy  f  A „n t-  2/\x/r,v  I  2.\s/n>,  !  2 A  ,^,v.  M) 

where  the  z\tJ  arc  the  component >  of  the  pioznspectroscopic  tensor  lor  the  defect.  The 
piczospcctroscopic  tensor  components  depend  on  the  symmetry  allowed  possible  defect 
orientations  in  the  crystal. 
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Thc  symmetry  of  the  'O'7  run  NaF  color  center  used  in  this  work  was  previously  examined 
with  polarized  spectroscopy  of  the  entire  /IM.  under  uniaxial  stress  '  and  pseudo-Stark 
splitting  of  PS  1 1  s  under  DC  electric  fields'1.  The  007  nm  NaF  color  center  has  (',  symmetry 
and  its  ZPF  transition  moment  a  is  perpendicular  to  the  [ll<>'!  oriented  defect  reflection 
planes  2°.  Defects  with  C,  symmetry  in  a  cubic  lattice  arc  referred  to  as  monoclinic  type  I 
defects,  for  which  there  are  twelve  equivalent  hut  distinct  orientations  of  the  {100}  oriented 
defect  planes13.  The  presence  of  a  reflection  plane  of  symmetry  limits  the  number  of 
independent  A«j  to  4,  which  arc  labelled  A,  with  i  =  1.2,7  and  4.  fable  I  presents  the  energy 
shift  coefficients  in  Hq.  4  for  each  of  the  twelve  defect  orientations.  By  examining  fable  1 
and  Fq.  4,  one  can  see  that  each  defect  orientation  has  a  dill'cicnt  absorption  frequency  shift 
for  an  arbitrary  stress.  With  no  applied  stress  or  strain,  defects  with  the  same  transition 
energy  arc  orientationallv  degenerate.  The  PS  1 1  splitting  results  when  phonon  stress  lifts  the 
orientational  degeneracy. 

To  use  these  frequency  shifts  as  a  detector  of  stress  or  strain,  we  must  also  know  the 
amount  that  each  of  the  possible  defect  orientations  contributes  to  a  PS! I,  i.e.  we  need  the 
size  of  the  components  in  the  splitting  pattern.  These  arc  not  the  same  sj/c  as  the 
components  in  a  linear  absorption  splitting  pattern.  Defects  with  transition  moments  in  a 
direction  /c  arc  hole  burned  at  a  rate  propoitional  to  la,*  pi  .  where  p  is  the  laser 
polarization.  I:or  shallow  holes,  the  contribution  to  the  hole  depth  is  proportional  to  the 
product  of  the  rate  at  which  defects  of  a  particular  orientation  arc  hole  burned  and  their 
unburned  absorption  strength,  i.e.  defects  with  /#,  contribute  an  amount  f-  |/i,»p|4  to  the 
hole  depth.  In  the  experiment,  the  laser  is  linearly  polari/ed  in  the  x-z  plane  (sec  l  ig.  I) 
with  direction  cosines  a  and  y  from  the  x  and  /  directions,  fable  I  also  shows  the 
unnormalized  contribution  f  that  each  defect  orientation  provides  to  the  hole  depth  in  terms 
of  the  laser  polarization  direction  cosines.  With  the  values  in  fable  I.  any  splitting  pattern 
may  be  computed  for  a  persistent  spectral  hole  burned  in  the  NaF  («07  nm  color  center  /IM. 
for  an  arbitrary  stress  and  laser  polarization  in  the  x-z  plane. 


Before  using  the  complete  hole  splitting  pattern  to  calculate  the  phonon  signals,  we 
consider  the  signal  that  results  from  a  net  frequenev  shift  of  the  entire  hole.  A  shift  of  the 
entire  hole  displaces  the  FM  spectra  in  Fig.  3  along  the  frequency  axis.  The  highest  sensitivity 
is  obtained  at  the  center  of  the  dispersion-like  S-  signal  in  Fig.  3  (I'd,  since  this  is  the  position 
of  highest  slope  in  the  PS  1 1  FM  spectra.  Note  that  the  sign  of  the  signal  indicates  the  sign 
of  the  net  frequency  shift  of  the  hole.  On  the  other  hand,  when  the  stress  produces  a 
symmetric  splitting  of  the  hole,  no  large  signal  results  for  laser  frequency  fixed  at  the  center 
of  the  S2  spectrum,  lienee,  it  is  also  useful  to  detect  the  effects  of  stress  for  another  laser 
frequency  and  RF  phase,  as  shown  in  Fig.  3  (a).  With  the  laser  frequenev  fixed  at  the  S, 
negative  peak  in  Fig.  3  (a),  any  splitting  pattern  of  the  hole,  symmetric  or  asymmetric, 
decreases  the  hole  depth  which  is  detected  as  an  increase  in  signal.  I  hc  sign  of  the  signal 
in  S|  phase  is  always  positive  regardless  of  the  type  of  splitting.  In  other  words,  the  center 
of  S2  is  used  to  phase  sensitively  detect  net  frequency  shifts  due  to  asymmetric  splitting,  and 
the  peak  of  S,  is  used  to  detect  any  splitting. 


The  S2  signal  for  the  three  probe  laser  polarizations  used  in  the  experiment  may  now  be 


computed.  The  dispersion  for  a  Lorent/ian-shaped  shaped  hole  may  be  written 
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where  Sp  is  a  magnitude  factor,  I  is  the  lull  width  at  half  maximum  of  the  absorption  hole. 


and  the  hole  is  centered  at  zero  frequency.  The  S,  spectrum  is  proportional  to 


22.  2.1 


S2(r)  ~  (f)(v  -  VK,  )  f  </>(>•  4  )  -  2r/>(\’). 


((') 


The  approximation  assumed  in  Fq.  6  is  the  limit  of  small  modulation  index  and  shallow 
holes.  A  modulation  index  near  1  is  used  in  these  experiments,  (or  which  higher  order  terms 
in  t>^  should  be  important  in  the  entire  spectrum  S,(rf  However,  uc  will  only  be  concerned 
with  the  slope  of  the  S2(v)  spectrum  near  the  center  of  the  I’SII  and  thus  can  neglect  these 
higher  order  terms  with  large  spacing  from  the  center  ol  the  hole,  for  a  small  shill  of  the 
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entirc  hole  by  A,  a  linear  approximation  to  the  the  S>  phase  signal  with  laser  frequency  fixed 
at  the  burn  frequency  is 


where  r(vRr,  T.  <5P)  is  constant  for  fixed  vRj.,  P.  and  <V  Pq.  7  reasonably  predicts  that  the 
slope  of  the  S2  spectrum  increases  with  the  hole  si/e  fip,  increases  with  the  PM  sideband 
spacing  vRF,  and  decreases  with  the  hole  width  T.  When  the  orientational  degeneracy  of  the 
hole  is  lifted  by  stress,  the  S2  signal  is  proportional  to 


where  the  component  sizes  f  and  the  shifts  A,  arc  computed  from  the  information  in  Table 
1.  Note  that  through  Table  1,  Pq.  8  is  linear  in  the  small  stress  or  strain  tensor  components 
of  a  phonon  (Pq.  8  may  be  written  in  terms  of  the  strain  using  Pq.  3). 

In  the  experiments  presented  in  subsequent  sections,  three  dilTcrcnt  probe  laser 
polarizations  where  examined:  (1)  +  45  0  (<*  ~  V  ~  l/v  2  ),  (2)  -45  (<z  =  —  y  —  l/x/T),  and  (3) 

horizontal  or  0  °  (a  =  1,  y  =  0).  The  S2  signals  for  these  three  laser  polarizations  arc 
computed  from  Pq.  8  and  Table  1  to  be 

S2+45 — -j|j-  [(A,  +  dA,Wxx  +  (8A,  I  2Aj)nn  MA,  i (P\,)r7v  -  lf>A  (Ti) 

S^45~y--j^-[(A,  TOA^K.  -MSA,  f  2A?KV  l  (A,  I  ('Mi) 

»-AiKv  I  4-  <  A I  l  A  :)a/f  .  (Me) 
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The  quantities  in  Ficjs.  9  involve  one  off  diagonal  or  ‘  hear  stress  component,  and  linear 
combinations  of  the  diagonal  or  axial  stress  components. 

Taking  sums  and  differences  of  the  S j4'  and  S;1'  expressions  gives 


1  c  +  J'\  1  *  \ 

y(s2  -s2  I'-pj'W 


final 


y  (sr*5  +  S245)  ~  T'lV  C(A'  +  <^\2K<  +  f8AI  1-  T\,KV  +  (A,  +9A2K7],  (10b) 


Notice  that  taking  the  sum  and  difference  of  the  Si'”  and  S: ”  signals  separates  out  one 
ofT-diagonal  shear  stress  from  the  diagonal  or  axial  stresses,  lit  other  words,  through  suitable 
choices  of  laser  polarization  and  linear  combinations  of  small  signals,  it  is  possible  to  obtain 
individual  stress  tensor  components  for  phonons  propagating  inside  a  crystal.  By  converting 
the  measured  phonon  stress  tensor  components  to  strain  with  l:q.  3.  I:q.  2  can  then  be  used 
to  obtain  information  about  the  phonon  wave  vectors  and  polarizations. 

The  high  symmetry  of  a  cubic  crystal  simplifies  computation  of  the  strains  from  stress. 
Cubic  symmetry  limits  the  form  of  the  clastic  stillness  tensor  in  !:q.  3  such  that  the  diagonal 
stress-strain  relations  arc  decoupled  from  the  oil-diagonal  stress-strain  relations15.  Tor 
instance,  the  ofT-diagonal  stress  and  strain  arc  simply  related  by  r„  —  n7tl2(\,,  i  j.  and  the 
diagonal  stress  and  strain  mix  among  themselves. 

Using  Kaplyansii's  uniaxial  stress  calculations,  and  Baumann's  uniaxial  stress  data. 14 
an  estimate  for  the  piezo  spectroscopic  components  may  be  obtained  for  the  NaT  607  nm 
color  center.  The  A(  arc  presented  in  Table  2.  These  cocllicicnts  will  be  useful  in  gauging 
the  sensitivity  of  the  technique  used  here  in  measuiim:  plmnon  stress  and  strain. 
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IV.  Phonon  detection  sensitivity  and  probe  polarization  dependence 


A.  Experimental  results 

Figure  4  shows  the  phonon  timc-of-flight  signals  lor  three  probe  laser  polarizations 
detected  in  S,  and  S2  I'M  phase  with  the  probe  laser  frequency  fixed  at  the  arrows  in  Fig. 

3.  The  probe  laser  was  positioned  at  x  =  0.04  cm  from  the  ('\  surface  and  was  vertically 
centered  (z  =  0).  The  0.1  cm  diameter  YA(i  laser  spot  was  horizontally  centered  on  the  Cr 
face  (y  =  0),  and  displaced  to  z  =  0.0.5  cm.  (The  importance  of  the  Y/\(i  spot  displacement 
will  become  apparent  in  Sec.  V).  The  YAG  laser  firing  time  is  at  zero  on  the  time  axis  in 
Fig.  4.  The  S2  phase  signals  are  at  0  V  and  the  S,  phase  signals  arc  near  -0.75  V.  The 
separation  between  the  S,  and  S2  signals  near  t  =  o  is  the  S,  I'M  spectrum  negative  peak 
size,  and  approximately  the  size  of  the  S2  I  M  spectrum  peak  value.  For  1500  averages  in 
a  20  iVtffz  bandwidth,  the  noise  level  is  about  5  %  of  the  peak  value  of  the  S,  or  S2  spectra, 
which  corresponds  to  a  change  in  transmission  of  about  0. 1  *V  The  signals  shown  in  Fig. 

4  clearly  represent  shifting  and  splitting  of  the  persistent  spectral  hole,  since  no  signals  arc 
obtained  if  the  same  experiments  arc  performed  with  no  persistent  spectral  hole  present. 

Signals  larger  than  the  noise  only  appear  more  than  1.5  // s  after  the  YAG  laser  fires. 
The  S,  signals  arc  relatively  simple  compared  to  the  S.  signals  --  each  of  the  S,  signals 
contains  two  large  peaks  with  some  structure.  The  S,  signals  show  no  qualitative  dependence 
on  probe  laser  polarization.  By  comparison,  the  S,  signals  show  much  more  complicated 
structure.  The  largest  S2  signals  correspond  to  the  largest  S,  signals,  but  there  arc  many 
more  Sa  signals  above  the  noise  in  f  ig.  4.  Some  signals  are  expected  to  be  smaller  in  S, 
phase,  since  the  PS  1 1  peak  size  is  less  sensitive  to  small  frequency  shifts  than  the  high  slope 
of  S2  (sec  Fig.  3),  For  purposes  of  later  dismission,  each  of  the  various  signals  has  been 
labelled  I.,,  I,2,  I.,,  T,.  T2.  Id.  and  Id  in  Fig.  4. 


The  features  !,,,  l,2,  l.j  and  l.j,  which  show  an  arrival  time  correspondence  in  the  S,  and 
S2  signals  can  be  compared  for  an  understanding  of  the  splitting  behavior  of' the  I’SII.  As 
described  in  See.  III.  the  S2  signals  only  respond  to  asymmetric  splitting  (with  a  net  shift) 
and  S,  records  any  splitting,  however  less  sensitively  than  S>.  The  I.i  pulse  shows  a  net 
upwards  frequency  shift  of  the  hole  in  any  polarization.  Comparing  the  S,  and  S2  signals 
as  I.,  crosses  over  into  I2.  one  observes  the  PSII  to  partially  recover  its  unsplit  condition, 
and  then  split  again  with  a  net  negative  frequency  shift.  The  I  ,  pulse  is  very  narrow  (20 
MHz  bandwidth  limited),  but  the  l.2  tail  following  the  I  ,  pulse  lasts  for  about  a  microsecond. 
Note  that  the  l.j  and  l.j  signals  mimic  the  I.,  and  1.2  signals,  but  with  an  overall  change  of 
sign,  and  a  shorter  l.j  tail. 

Unlike  the  I.  signals,  the  T  signals  show  a  qualitative  dependence  on  the  laser 
polarization.  Note  that  the  f  signals  in  I  ig.  4  (a)  and  (b)  each  consist  of' a  frequency  swing 
of  one  sign  followed  by  a  frequency  swing  of  the  opposite  sign.  This  shape  shall  subsequently 
be  referred  to  as  the  "derivative  shape"  of  the  I  signals.  The  derivative-shaped  features 
labelled  Ti  and  T2  change  sign  for  probe  laser  polarization  changing  from  -M5  0  to  -45  °  in 
Pig.  4  (a)  and  (c).  Note  that  the  T,  and  T?  signals  arc  much  smaller  in  Pig.  4  (b)  for 
horizontal  probe  laser  polarization. 

One  perplexing  aspect  of  the  S2  signal  is  its  richness.  A  simple  view  of  pulsed  phonon 
propagation  is  that  for  an  arbitrary  displacement  within  the  solid,  one  longitudinal  and  two 
transverse  modes  arc  generated,  with  the  transverse  modes  propagating  with  roughly  half  the 
longitudinal  speed  in  NaP.  The  multitude  of  S,  arrival  times  observed  cannot  all  be 
reconciled  with  phonons  travelling  in  nearly  a  •  1 00  •  crystal  direction.  An  identification 
of  the  phonons  giving  rise  to  these  signals  will  be  aided  bv  the  laser  beam  position 
dependences  in  the  Section  V.  In  the  next  subsection,  the  types  of  stress  and  strain  that  give 
rise  to  the  signals  in  Pig.  4  arc  examined  using  the  results  of  Section  III. 
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B.  Analysis:  stress  and  strain  detection  sensitivity 

V  * 

To  determine  the  I'M  signal  sensitivity  to  phonon  stress  and  strain,  an  individual  stress 
component  will  be  separated  out  as  described  by  Fqs.  10.  Figure  5  shows  linear  sums  and 
differences  of  the  Sj45  and  S24?  traces  taken  from  Fig.  4.  (a)  and  (b).  Fig.  5  (a)  is  proportional 
only  to  the  time  variation  of  the  n,x  stress  or  strain.  Figure  5  (b)  is  proportional  »o  a 
linear  combination  of  diagonal  stress  or  strain  r.{l  components.  The  S?  signal  from  Fig.  4 
(b)  is  presented  in  Fig.  5  (c)  for  comparison.  Comparison  of  Fig.  5  (h)  and  (c)  shows  they 
are  qualitatively  similar,  which  is  consistent  with  the  idea  that  each  is  a  linear  sum  of  n,*  (Fqs. 
10b  and  9c).  Before  computing  the  sensitivity  to  stress  and  strain,  a  few  initial  comments 
about  the  varieties  of  stress  and  strain  giving  rise  to  the  I,  and  I  signals  shall  be  made,  which 
will  be  used  in  phonon  identification  in  later  sections. 

The  largest  features  in  the  r.„  trace  in  Fig.  5  (a)  are  the  T,  and  f2  derivative  shapes. 
Notice  that  the  T,  and  T2  derivative  shapes  are  missing  from  l  ig.  5  (h).  Instead,  the  T,  and 
T:  signals  in  Fig.  5  (b)  each  consist  of  broad,  positive  frequency  swings.  In  other  words,  the 
derivative  shape  in  the  T(  and  T2  signals  results  only  from  r.,x  shear  strain.  Also  note  that 
the  T,  and  T2  derivative  shapes  arc  about  the  same  si/e  in  Fig.  5  (a),  but  in  Fig.  5  (b)  the 
T,  positive  peak  is  larger  than  the  T2  peak.  Thus,  the  T,  and  I :  signals  have  similar  r.,x  shear 
strain,  but  dissimilar  sums  of  r.„  axial  strain  components. 

In  contrast  with  the  T  signals  in  Fig.  5  (a),  the  largest  signals  in  Fig.  5  (b)  and  (c)  arc 
the  F  features.  Notice  that  the  F  signals  arc  largely  missing  from  Fig.  5  (a),  lienee,  the  F 
signals  result  only  from  axial  stress  n„  or  strain  r.„  components. 

The  sensitivity  of  the  technique  to  an  individual  stress  or  strain  component  is  computed 
as  follows.  With  F500  averages  and  a  20  Mil/  bandwidth,  the  S2  signals  have  a  noise  level 
which  is  5  %  of  the  peak  of  the  S2  spectrum,  or  about  28  mV  in  Fig.  4.  The  slope  at  the 
center  of  the  S2  spectra  for  Fig.  4  is  found  using  the  slope  in  Fig.  2  (b)  and  the  S,  to  S2  signal 
spacing  in  Fig.  4  to  be  about  14  mV/ Mil/.  (Computing  the  prefactor  F/F  in  Fq.  7  gives 
about  10  mV/MIF/  if  one  uses  0.7.2  V  as  a  conservative  guess  Ibr  dr  and  F  —  2.20  Mil/.  Wc 
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will  use  F/r  =  l4mV/MIIz  below.)  The  noise  level  in  the  combination  of  signals 
(Sj45  -  Sj't5)/2  in  Fig.  5  (a)  is  38/V 2  =  27  mV.  Hence,  a  signal  to  noise  ratio  S/N  =  I  in 
(Si45  —  Si45)/2  corresponds  to  an  average  frequency  shift  of  the  PS! I  spectrum  bv  27/14  =  2 
MHz,  which  coincidentally  is  the  jitter  line  width  of  the  laser.  Htili/ing  f-'/F  =■  l4tnV/MHz 
to  convert  the  proportionality  in  Fqs.  9  and  10  to  an  equality,  the  n,x  producing  S/N  =  I 
in  Fig.  5  (a)  is  a,x  =  220  N/mJ  =  0.002  atmospheres.  The  corresponding  shear  strain  is 
c7,  =  ozxl 2Cm  =  4  x  10  9. 

One  should  keep  in  mind  that  the  values  computed  above  arc  averages  over  the  probe 
volume,  which  the  phonons  may  not  uniformly  irradiate.  In  general,  crystalline  elastic 
anisotropy  causes  the  phonon  group  velocity  and  the  phase  velocities  to  have  different 
directions.  A  large  solid  angle  of  phase  velocity  is  often  channeled  into  a  small  solid  angle 
of  group  velocities  resulting  in  an  energy  (lux  enhancement  in  certain  crystal  directions.  This 
phenomenon  has  been  termed  phonon  focusing16.  Fnergy  focusing  factors  for  special  crystal 
directions  have  been  computed  by  Maris,  16\vhich  give  values  between  3  and  0.4  for  NaF  at 
low  temperatures.  Elastic  wave  energy  is  quadratic  in  the  stress  and  strain,  lienee,  for  S/N 
=  I  in  Fig.  5  (a)  the  corresponding  shear  stress  and  strain  fields  inside  the  probe  volume 
may  have  peak  magnitudes  larger  than  220  N/m’  and  4  x  |0  4  by  a  factor  of  v'3  . 


These  results  may  be  compared  and  contrasted  with  the  earlier  phonon  detection  experiments 
of  Tcngfcllncr,  ct  al.14.  The  acoustic  echoes  detected  in  that  work  by  observing  the  changes 
in  transmission  at  the  center  of  the  spectral  hole  were  not  sensitive  to  the  phase  of  the 
ultrasonic  pulse.  The  sensitivity  of  our  I  MS  approach  is  sufficient  to  exceed  the  strain 
sensitivity  of  the  earlier  work  with  an  increased  time  resolution  (from  5  /is  to  50  ns),  and 
allow  for  phase  sensitive  detection  of  phonons  up  to  20  MHz.  While  this  simple  comparison 
has  not  precisely  included  the  difference  in  detection  bandwidth  for  the  two  experiments,  it 
is  clear  that  at  least  the  phase-sensitivity  of  the  IMS  approach  is  unique. 
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V.  Laser  position  dependence  and  phonon  identification 


A.  Experimental  results 

To  further  investigate  the  origin  of  the  various  signals  in  Figs.  4  and  5,  the  pump  and 
probe  beam  positions  arc  varied.  Figure  6  shows  YACi  spot  position  dependence  with  a  0. 1 
cm  diameter  YACi  spot  horizontally  centered  (y  =  0),  but  with  different  vertical  positions 
(different  z).  The  YACi  spot  z  positions  arc  indicated  on  the  side  of  Fig.  6.  I'hc  I'M  probe 
beam  is  vertically  centered  (z.  =  0),  and  set  near  the  middle  of  the  crystal  at  x  =  1.07  cm. 
t  or  the  traces  in  Fig.  6,  a  single  persistent  spectral  hole  was  burned  and  probed  with  +  45° 
laser  polarization.  By  examining  Fig.  6.  one  observes  that  the  I.  signals  depend  little  on  the 
YACi  spot  vertical  position.  By  contrast,  the  derivative  shaped  T  signals  show  a  strong 
dependence  on  the  YACi  spot  position. 

As  the  YACi  spot  moves  through  the  probe  laser  plane  z  =  0.  the  derivative  shaped 
portion  of  the  T  signals  become  narrower,  disappear  for  z  =  0.  and  reappear  with  the 
opposite  sign  and  broaden  again.  When  the  YACi  spot  is  centered  in  the  plane  of  the  probe 
laser,  the  T  signals  that  remain  appear  to  be  the  positive  frequency  swings  found  in  Fig.  5 
(b).  In  other  words,  the  positive  and  negative  parts  of  the  n.r  derivative  shaped  f  signals 
overlap  in  time  and  cancel  when  the  lasers  arc  vertically  centered  in  the  crystal.  A  possibility 
for  the  change  in  sign  of  the  arx  derivative  shape  is  that  two  different  stress  pulses  containing 
n,x  of  opposite  sign  travel  over  different  paths  to  the  probe  beam.  With  the  YACi  laser 
centered,  these  path  lengths  balance  and  the  combined  n  x  stresses  cancel.  With  the  YAC? 
beam  vertically  off  center  the  paths  are  not  the  same  length  and  the  n,x  of  one  sign  arrives 
before  the  a,t  of  the  other  sign.  In  other  words,  these  paths  involve  some  vertical  (z) 
component  of  motion  so  that  one  path  becomes  shorter  and  the  other  longer  by  vertical 
YACi  displacement.  To  further  investigate  the  different  paths  that  the  phonons  follow,  the 
probe  beam  is  also  displaced  within  the  crystal. 


Figure  7  shows  the  pulsed  phonon  signal  dependence  on  the  probe  beam  horizontal 
position.  The  YAG  laser  spot  is  centered  on  the  crystal  end  lace.  I  he  probe  beam  is 
vertically  centered  (z  =  0),  and  displaced  horizontally  in  increments  of  about  0.1.5  cm.  The 
probe  beam  x  values  arc  shown  on  the  side  of  Fig.  7.  For  each  probe  laser  position,  a  new 
persistent  spectral  hole  is  burned  and  probed  with  t- 45°  probe  laser  polarization.  Hashed 
lines  arc  drawn  in  Fig.  7  to  mark  the  movement  ol'  the  various  signals  as  the  probe  beam 
moves. 

Considering  the  F  signals  for  the  moment,  one  sees  that  F,  and  1,2  arrive  later  and  FJ 
and  FJ  arrive  earlier  as  the  probe  laser  moves  away  from  the  Cr.  It  is  apparent  that  F,  and 
F2  arc  two  parts  of  a  single  signal,  and  similarly  lor  FJ  and  I  The  F,- F2  arrival  time  is 
consistent  with  longitudinal  acoustic  (FA)  phonons  traveling  in  a  [100]  crystal  direction  at 
(6.24  + 0.03)  x  10s  cm/s,  which  is  the  highest  acoustic  velocity  in  NaF.  As  the  probe  laser 
moves  away  from  the  Cr.  it  moves  closer  to  the  other  end  of  the  crystal.  The  decreasing 
FJ-F]  arrival  time  is  consistent  with  the  return  of  the  FA  [1 00]  phonons  to  the  probe  laser 
volume  after  a  reflection  from  the  far  crystal  surface.  In  longer  time  traces  (not  shown  here), 
as  many  as  10  reflected  FA  [100]  phonon  pulses  arc  observed,  with  the  time  interval  between 
every  other  pulse  corresponding  to  an  FA  [1 00]  phonon  round  trip  time  in  the  crystal.  As 
observed  in  Fug.  7,  the  reflected  FJ-F)  signal  is  opposite  in  sign  to  the  FrF2  signal.  After 
each  reflection  a  sign  change  is  observed,  so  that  FA  [100]  phonons  travelling  in  the  +x 
direction  have  the  Ft-F2  sign,  and  reflected  FA  phonons  travelling  in  the  -x  direction  have 
the  FI-Fj  sign.  This  180°  phase  change  is  expected  for  an  acoustic  wave  reflected  at  an 
interface  of  low  acoustic  impedance  to  high  acoustic  impedance  such  as  the  NaF-lIc  interface 
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The  movement  of  two  other  signals  can  clearly  be  followed  in  Fig.  7.  The  F,  and  f, 
signals  also  arrive  at  later  times  as  the  probe  laser  moves  away  from  the  Cr  surface,  but 
more  slowly  than  at  a  fixed  multiple  of  the  I  ,1  ,  I  A  f  100]  arrival  time.  Data  similarly  taken 
with  the  YAG  spot  off  center  show  that  the  >  T,  derivative  shape  (see  Fig.  6)  follows  the 
movement  of  the  ru  I]  peak  in  Fig.  7.  It  is  difficult  to  follow  the  movement  of  the  l\  signal. 
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as  the  region  after  the  T,  signal  becomes  compiicntcil  lor  some  probe  laser  positions.  A 
possible  reason  for  this  will  be  explored  in  the  next  section. 

The  data  in  Fig.  7  arc  presented  in  a  different  wav  in  Fig.  ,X  to  aid  further  in  signal 
identification.  In  Fig.  8,  the  time  axis  for  each  trace  in  Fig.  7  has  been  normalized  to  the 
F,  pulse  arrival  time,  which  transforms  the  time  axis  to  a  relative  velocitv  axis  for  [100] 
propagating  phonons.  In  particular.  [100]  propagation  transverse  acoustic  (  FA)  phonons 
which  travel  at  approximately  1/2  the  FA  [ion]  speed  should  produce  a  signal  near  2  on  the 
relative  time  axis  in  Fig.  8.  Remarkably,  an  examination  of  Fig.  8  shows  that  there  .arc  no 
additional  signals  that  line  up,  i.c.  the  I  A  [loo|  phonons  appear  to  be  missing.  Signals  that 
do  not  line  up  when  the  time  axis  is  converted  to  a  relative  speed  axis  as  in  l  ig.  8  must  either 
(I)  result  from  delayed  emission  from  the  surface,  or  (2)  do  not  travel  in  [100]  crystal 
directions,  or  (3)  do  not  travel  ballistic-ally  (at  constant  speed).  I  he  present  data  cannot  be 
used  to  rule  out  possibilities  (I)  or  (3),  but  condition  (2)  is  the  most  reasonable  possibility 
consistent  with  the  results  of  Figs.  6  and  7. 

In  Fig.  6,  it  was  found  that  different  paths  arc  followed  by  the  positive  and  negative  parts 
of  the  r.7X  pulses  in  the  Ti  and  T2  signals.  Also,  these  paths  involve  some  vertical  component 
of  motion  so  that  one  path  becomes  shorter  and  the  other  longer  by  vertical  YA(i 
displacement.  This  is  consistent  with  one  possibility  from  I  \v.  S  that  the  F,  and  l2  signals 
arc  not  from  TA  phonons  propagating  in  nearly  the  [  ion]  crystal  direction.  In  order  to  have 
both  a  vertical  component  of  motion  and  to  reach  the  probe  volume,  these  phonons  must 
reflect  off  the  upper  and  lower  side  walls  of  the  crystal.  One  may  then  propose  that  the 
Ti,  1 1  and  T2  signals  all  result  from  phonons  propagating  in  directions  widely  dilfcrent  from 
the  [100]  crystal  direction.  In  the  next  section,  an  identification  of  the  I  I and  I  ?  signals 
as  side  wall  reflected  phonons  is  made. 
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B.  Analysis:  Phonon  identification 

Ihrcc  different  concepts  about  phonon  generation  and  propagation  in  crvstals  are  used 
in  this  section  to  explain  the  various  phonon  stress-strain  signals  observed  in  Nab.  The  main 
results  are  reviewed  below,  and  then  discussed  in  erms  of  t lie  phonon  hot  spot  at  the  surface, 
phonon  focusing  in  crystals,  and  phonon  mode  conversion  upon  reflection.  The  features 
observed  by  measuring  the  stress  and  strain  of  pulsed  phonons  inside  of  a  N'al  crystal  arc 
summarized  as  follows.  I  lie  I.A  [  1 00 J  propagating  phonons  are  the  first  phonons  to  arrive 
and  consist  of  a  comprcssional  pulse  followed  hv  a  dilation.il  tail  (I.,  --  f.:).  TA  [100] 
phonons  arc  expected  to  arrive  at  nearly  twice  the  I.A  [loo]  phonon  arrival  time,  but  arc 
not  observed  by  our  choice  of  geometrv  and  polarizations.  Instead  of  observing  the  TA 
[100]  phonons,  many  other  phonon  signals  ate  observed  ( I ...  I',,  and  I',),  and  their  behavior 
is  consistent  with  side  wall  reflected  phonons  travelling  in  directions  widely  different  from  the 
[100]  crystal  direction.  The  shear  strain  and  r„  axial  strain  portions  of  the  T,  and  Tj 
signals  have  been  separated,  as  in  I  ig.  5.  I  he  f.  and  I  -  derivative  shapes  result  from  r.t 
shear  strains  and  arc  about  the  same  in  si/e  (Tig.  5  (a)).  I  he  r,,  strains  reflected  from  the 
upper  crystal  face  arc  opposite  in  sign  to  those  reflected  from  the  lower  crystal  face  and  the 
difference  in  arrival  times  for  an  off  center  YA(i  spot  produces  the  positive  and  negative 
going  derivative  shape.  I  he  sum  of  axial  t„  strain  components  for  |  ,  and  l\  appears  to  be 
comprcssional,  and  the  T,  peak  is  larger  titan  that  for  I.-  (fig.  5  (hi).  The  small  I  ,  negative 
frequency  shift  of  the  hole  appears  to  involve  a  dilation al  sum  of  r„  axial  strains. 

In  order  to  discuss  phonon  propagation  within  a  crystal  and  mode  conversion  upon 
reflection,  it  is  necessary  to  introduce  t lie  slowness  and  group  vcfocitv  surfaces,  figure  0 
shows  the  slowness  and  group  vcfocitv  mii fares  computed  for  \a|  at  d  K  in  a  (001)  crystal 
plane.  These  have  been  computed  using  published  Immulac  and  -1  K  clastic  constants  ”, 

I  he  slowness  vector  is  in  the  same  direction  as  the  phase  velocity  or  k  but  has  ,i  magnitude 
that  is  the  reciprocal  of  the  phase  velocity  I  here  ate  three  mu  h  sutlaces  (I  ig.  0(a)).  |  hc 

outermost  surface  with  circular  cross  section  in  the  (0<t|)  plane  is  the  pure  shear  or  pure 
transverse  (FT)  surface,  and  has  polarization  imintal  to  the  Htt'l)  plane.  Ihc  surface 


degenerate  with  the  PI  surface  in  100  '  crysial  diicctions  is  the  quasitransversc  (QT) 


surface,  and  has  polarization  in  the  plane,  hut  not  noccssaiih  normal  to  the  slowness  or  k 
direction,  flic  innermost  surface  is  the  quasilongiUulinal  (Ol  )  surface,  with  polarization  in 
the  plane  and  nearly  parallel  to  k.  A  point  on  the  surface  specifics  the  direction  of  the 
slowness  or  the  wave  vector  k.  As  shown  in  I  ig.  0  (at.  the  normal  to  the  slowness  surface 
specifics  the  direction  of  the  corresponding  group  velocity  V,..  Vote  that  the  group  velocity 
and  phase  velocity  arc  not  necessarily  in  the  same  direction.  Most  of  the  normals  to  the 
QT  surface  point  in  directions  close  to  the  --  I  10  >  directions.  In  other  words,  large  solid 
angles  of  k  for  QT  phonons  are  channeled  imo  small  solid  angles  of  group  velocities  near 
the  •  I  10.'*  crystal  directions.  Phis  channeling  of  the  energy  into  a  small  solid  angle  is 
known  as  phonon  focusing  ".  The  concave  shape  of  the  QT  surface  produces  the  cusps  in 


the  group  velocity  surface  in  fig  ')  (h).  and  these  arc  a  characteristic  signature  of  phonon 


focusing 
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Inspection  of  the  QT  surface  in  l  ig.  0  (a)  shows  that  there  is  some  enhancement  of  the 
QT  phonons  not  in  the  <  II0>  directions,  hut  rather  in  the  ■  100  *  directions,  focusing 
factors  expressing  the  amount  of  enhancement  in  high  symmetry  directions  have  been 
calculated16.  Using  4  K  clastic  constants  1  to  compute  the  enhancement  factors  for  Naf, 
one  finds  QT  phonons  arc  enhanced  hv  a  factor  of  TO  near  the  <  |O0>  directions,  and  QT 
phonons  arc  enhanced  by  a  factor  of  2.3  near  the  •  110  crystal  directions.  I  nhanccmcnt 
of  the  the  TA  phonons  in  the  <"  1 00  >  crystal  directions  is  consistent  with  the  strong  l.A 
[|00]  phonon  signal  labelled  Tt  —  T».  An  explanation  for  the  absence  of  [100]  propagating 
1  A  phonon  signals,  which  would  have  been  measured  by  their  r„  trains,  is  that  the  TA 
phonons  arc  dcfocuscd  from  the  [100]  direction.  Iloweset.  I  A  phonons  are  focused  in 
"■  I  10  >  directions,  which  gives  rise  to  the  \aiicty  of  side  wall  lellected  phonons  as  explained 
below. 


Phonon  focusing  considerations  are  consistent  with  the  observation  of  a  strong  I  A 
[100]  phonon  signal,  no  detectable  I  A  [loo]  phonon,  and  most  other  signals  appearing  as 
side  wall  reflections  if  one  takes  into  account  mode  (  (inversion  upon  reflection.  As  is 


well-known,  reflection  and  transmission  of  a  single  acoustic  mode  at  an  interface  generally 
involves  acoustic  mode  conversion  11  W.  The  I  I  ,,  and  I  .  signals  mav  then  he  explained 
as  follows:  I.A  phonons  incident  on  t he  side  walls  (travelling  in  approximately  the  110  > 
crystal  directions)  arc  specularly  reflected  as  !  A  phonons  and  partially  converted  to  TA 
phonons.  TA  phonons  incident  on  the  side  walls  arc  rellccted  specularly  as  TA  phonons  and 
partially  converted  to  I.A  phonons.  The  arrival  times  of  the  I  i  and  I  j  signals  arc  consistent 
with  LA  to  I.A  and  I  A  to  I A  side  wall  reflected  phonons,  respectively.  Strong  focusing  of 
the  I'A  phonons  in  -•  1 10  >  directions  is  responsible  for  the  specularly  reflected  TA  to  TA 
l2  signal  appearing  strongest  when  the  probe  laser  beam  is  centered  horizontally  in  the  x 
direction,  as  in  figs.  4  and  5.  I  he  complexity  in  the  signals  that  arivc  after  passage  of  the 
T,  signal  in  Figs.  7  and  S  can  probably  be  attributed  to  a  weaker  distribution  of  phonon 
modes  that  arc  non-spccularlv  reflected  upon  mode  conversion.  The  arrival  time  of  the 
intermediate  T,  signal  is  consistent  with  both  I  A  to  TA  conversion  and  TA  to  I.A  conversion 
upon  reflection. 

I  he  sign  and  magnitude  of  the  strains  belonging  to  the  observed  TA  phonons  in  T,  and 
I2  arc  also  consistent  with  the  side  wall  reflection  and  mode  conversion  identification.  T, 
and  Tj  signals  both  consist  of  shear  strains  that  travel  to  the  probe  volume  after  reflection 
from  upper  and  lower  side  walls  with  opposite  signs  for  >: ,  (I  ig.  5  (a)).  By  examining  l:q. 

2.  one  sees  that  TA  phonons  that  have  r.„  must  be  <)  I  types  of  phonons,  with  u  and  k  both 
containing  components  in  the  7-x  plane.  By  performing  the  coordinate  transformation 
7. —  7  one  obtains  r.  =  —  r.„.  i.c.,  by  turning  the  crystal  ISO  "  about  the  x  axis  the  sign 
of  the  7.x  shear  strain  changes.  This  is  consistent  with  the  r  ,  shear  component  of  the  strain 
pulse  reflected  from  the  upper  surface  being  opposite  in  sign  to  the  pulse  reflected  from  the 
lower  surface,  fn  f  ig.  5  (h),  the  larger  sum  of  for  I ,  than  T,  result  from  T,  consisting 
of  both  mode  converted  LA  and  I  A  phonons,  whereas  L  results  only  from  TA  to  I  A 
reflections  without  mode  conversion. 

The  final  result  to  be  discussed  is  the  appearance  of  t lie  dilation, al  tail  following  the 
compressions!  [100]  propagating  I.A  phonon  packet.  Heat  pulses  gcnctaied  at  surfaces  with 


pulsed  lasers  have  been  known  to  produce  a  phonon  hot  epot  consisting  of  a  broad 
spectrum  of  phonon  excitations'*".  The  experimental  manifestation  of  the  hot  spot  is  a 
tail  appearing  after  the  initial  phonon  pulse  as  detected  with  a  phase  insensitive  bolometer, 
tunnel  junction,  or  phonon  induced  fluorescence.  The  explanation  for  the  tail  is  that  the 
clastic  scattering  rate  for  phonons  is  exponential  in  the  phonon  frequency,  which  tends  to 
trap  high  frequency  phonons  near  the  surface,  l  he  high  frequency  phonons  decay 
anharmonically  to  lower  frequency  phonons,  which  then  escape  from  the  surface  region  at 
a  later  time.  With  our  phase-sensitive  technique,  we  observe  directly  that  the  initial  pulse 
is  compressions!  in  nature,  and  the  hot  spot  tail  is  dilations!.  This  new  observation  within 
a  single  crystal  at  low  temperatures  is  consistent  with  the  hot  spot  interpretation  if  thermal 
expansion  and  pressure  relief  arc  associates!  with  the  hot  spot,  as  has  been  surmised  from 

,  .11.  <14 

room  temperature  experiments  on  metals. 

VI.  Conclusion 

f’hase-scnsitive  optical  detection  of  pulscd-lascr-gcncratcd  stress-strain  waves  inside  of  a 
solid  has  been  demonstrated  using  frequency-modulation  spectroscopy  (FMS)  of  a  narrow' 
persistent  spectral  hole  ( PS  1 1 )  in  a  Nab  color  center  /cm  phonon  line  (XI’I).  lhe  inherent 
sensitivity  of  the  FMS  technique  and  the  spectral  sharpness  of  I’Sfl's  allowed  detection  of 
strains  as  small  as  4  x  10  9  in  a  20  Mil/  bandwidth,  focusing  the  probe  beam  diameter 
inside  the  NaF  crystal  to  less  than  one-half  the  minimum  acoustic  wavelength  expected  lor 
20  Mil/  phonons  in  NaF  allows  the  probe  beam  to  follow  the  phase  or  sign  of  phonons 
up  to  the  bandwidth  limit. 

High  sensitivity  to  strain  led  to  a  rich  vauetv  of  observable  phonon  signals  after  pulsed 
laser  generation.  Identification  of  the  various  types  of  phonon  wave  packets  was 
accomplished  by  selecting  different  RF  phases  of  detection,  by  burning  and  probing  I’SIFs 
with  different  laser  polarizations,  and  by  moving  the  phonon  generation  and  probe  laser 
beams  relative  to  the  crystal.  A  propagating  stress-strain  wave  was  shown  to  shift  the 
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frcqucncics  of  differently  oriented  color  center';  by  unequal  amount';.  By  varying  t he  probe 
laser  polarization,  which  varies  the  contribution  from  different  color  centers  to  the  detected 
signal,  it  was  possible  to  distinguish  between  propagating  axial  and  shear  strain  components. 
The  directions  of  propagation  of  the  phonon  packets  were  identified  by  displacing  the  phonon 
generation  and  probe  laser  beams.  A  combined  analysis  of  the  travelling  strain  types,  signs, 
and  propagation  directions  led  to  an  identification  of  the  various  phonon  packet 
polarizations. 

A  direct  extension  to  this  work  could  include  observations  of  surface  generated  phonon 
packet  propagation  in  other  crystal  directions  so  that  phonon  focusing  occurs  along  the  direct 
path  to  the  probe  region.  For  example,  with  a  [110]  oriented  crystal.  TA  phonon 
propagation  could  be  studied  without  side  wall  reflection  and  accompanying  mode 
conversion.  Since  I’SII  burning  is  found  to  occur  in  the  defect  absorption  spectra  of  many 
different  types  of  solids15,  the  phase-sensitive  phonon  detection  technique  demonstrated  here 
should  also  be  applicable  to  the  study  of  diffusive  phonon  propagation  n  glasses  and 
polymers,  as  well  as  to  ballistic  phonon  propagation  in  other  crystals 
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Figure  Captions 


Figure  1.  Sample  configuration  for  laser  generation  and  phase-sensitive  timc-of-flight 
detection  of  phonons.  A  thin  film  of  Cr  is  evaporated  on  one  face  of  a  single  crystal  of 
NaF.  The  crystal  contains  color  centers,  has  polished  !  lot))  laces,  and  is  immersed  in 
superfluid  He  at  1.6  K.  A  pulsed  Nd:YAG  laser  beam  is  used  to  generate  heat  pulses  in  the 
Cr  film.  A  probe  laser  beam  is  used  to  burn  a  narrow  persistent  spectral  hole  in  a  color 
center  zero  phonon  line.  Subsequent  frequency-modulation  of  the  probe  laser  beam  allows 
broadband  detection  of  the  modulation  of  the  hole  shape  by  the  phonon  pulses. 

Figure  2.  Experimental  apparatus  for  phase-sensitive  timc-of-flight  detection  of  phonons. 
The  subsystems  of  the  apparatus  arc  a  frequency-modulated  dye  laser  probe  beam  (DT).  a 
pulsed  Nd:YAG  laser  beam  (YAG)  used  to  generate  phonons  in  a  sample  (S)  held  at  1.6  K 
in  an  optical  immersion  cryostat,  and  standard  high  speed  laser  I'M  detection  electronics. 

A  detailed  description  of  the  apparatus  is  given  in  the  text. 

Figure  3.  Persistent  spectral  hole  frequency-modulation  spectra.  A  persistent  spectral  hole 
(PSH)  is  burned  into  the  NaF  607  nm  color  center  zero  phonon  line  at  1.6  K  while  the  laser 
frequency  modulation  is  off.  FM  spectra  of  the  PSII  are  measured  by  tuning  the  laser 
frequency  through  2  GHz  about  the  burn  frequency  while  the  laser  frequency  is  modulated 
at  2.50  MHz.  The  FM  spectra  arc  shown  in  (a)  S,  or  absorption  phase,  and  (b)  S2  or 
dispersion  phase.  The  arrows  point  to  the  frequencies  at  which  the  laser  is  fixed  in  order 
to  detect  the  modulation  of  the  hole  shape  bv  phonons. 

Figure  4.  Phonon  timc-of-flight  signals  for  three  probe  laser  polarizations.  The  probe  laser 
polarization  is  fixed  at  (a)  +45  °  from  horizontal,  (b)  horizontal,  or  (c)  -45  °  from  horizontal. 
Persistent  spectral  holes  (PSH)  arc  burned  and  probed  in  the  same  polarization.  For  a  given 
probe  polarization,  the  time  dependence  o(  the  PSII  is  monitored  at  the  negative  peak  of  the 
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Si  phase  [rM  spectrum  or  the  center  of  the  S»  phase  1  \ I  spectrum,  as  shown  in  l  ip.  3.  I'ach 
trace  is  an  overlay  of  two  traces,  each  of  which  is  an  average  over  1500  Nd:YAG  laser  shots. 

Figure  5.  Combinations  of  phonon  timc-of-flight  signals  lor  different  probe  laser 
polarizations.  Traces  (a)  and  (b)  arc  combinations  of  the  S7  phase  signals  shown  in  Fig.  4. 
(a)  The  difference  of  S2  phase  signals  for  probe  laser  polarizations  at  +45  "  and  -45  °  (i.c., 
the  difTcrencc  of  S2  traces  in  Fig.  4  (a)  and  4  (b)).  (b)  The  sum  of  +45"  and  -45°  polarization 
Sj  traces  from  Fig.  4  (a)  and  (b).  (c)  is  the  S2  phase  signal  measured  with  horizontal  (0") 
probe  polarization  for  comparison  with  the  combinations  in  (a)  and  (b). 

Figure  6.  NdtYAG  laser  spot  position  dependence  of  the  S2  phase  phonon  timc-of-flight 
signals.  For  each  trace,  the  0.1  cm  diameter  NdtYAG  laser  spot  is  horizontally  centered  on 
the  crystal  end  face  (y  =  0),  and  set  at  a  different  vertical  position  (z),  as  indicated  on  the 
side  of  the  figure.  The  probe  laser  position  is  lived  in  the  middle  of  the  crystal,  and  the  probe 
polarization  is  set  at  +4.5  °  from  horizontal.  Fach  trace  is  an  overlay  of  at  least  two  traces 
(five  traces  for  z  =  0,  three  taken  before  and  two  after  the  other  traces),  each  of  which  is  an 
average  of  over  1500  Nd:Y’AG  laser  shots. 

Figure  7.  Phonon  timc-of-flight  signal  dependence  on  probe  laser  position.  For  each  trace, 
the  probe  laser  is  positioned  inside  the  crystal  at  a  new  distance  from  the  C'r  surface,  and  a 
new  persistent  spectral  hole  is  burned  with  +45°  probe  polarization,  flic  probe  distance 
from  the  Cr  surface  is  indicted  on  the  side  of  the  figure.  I  he  phonon  timc-of-flight  signals 
arc  measured  in  S2  phase.  The  YAG  laser  firing  time  is  at  zero  on  the  time  axis.  The  dashed 
lines  in  this  figure  trace  the  movement  of  the  phonon  signals  as  the  probe  position  changes. 

Figure  8.  Phonon  timc-of-flight  signal  dependence  on  probe  laser  position  with  normalized 
time  axis.  The  data  in  Fig.  7  arc  presented  in  this  Fie.  with  the  time  axis  for  each  trace 
divided  by  the  F,  signal  arrival  time. 
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Figure  9.  Slowness  and  group  velocity  surfaces  for  \'al  at  4  K.  (a)  flic  acoustic  wave 
slowness  surfaces  arc  shown  in  the  (00 1 )  crystal  plane.  The  fiducial  marks  in  (a)  .arc  spaced 
1  x  10 ■*  s/ctn  apart.  The  normal  to  a  slowness  surface  gives  the  direction  of  the  group 
velocity  V,.  (b)  The  corresponding  group  velocity  surfaces  are  shown  in  the  (00 1)  plane. 
The  fiducial  marks  in  (h)  arc  spaced  2x  10*  cm/s  apart. 


Table  1.  Persistent  Spectral  Hole  Stress  Shift  Coefficients  for  Monoclinic 

(Type  I)  Defects  in  Cubic  Crystals 
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